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1. INTRODUCTION 

Vessels are sailing in bad weather conditions, due to effect of the external disturbances such as 
wave, wind and current can produce large amplitude roll motion, this is a dangerous phenomenon, it affects 
to the comfortable of the crew, cargo and ship, this phenomenon also affects to the vessels economic 
efficiency, therefore, research to reduce roll motion is a very important task. To perform that issue, there are 
some devices which have been researched to reduce roll motion [1]-[4] and can be named as bilge keels, 
antirolling tanks. But, due to the limitation of those devices, the active fin stabilizer has been researched with 
high efficiency in the roll reduction. Due to the superior of roll reduction ability, active fin has been 
researched very early and became a hotspot in the ship control filed [5]. Numerous researches on ship roll 
reduction by using active fin have been researched since 1940s. For example, a robust fin roll reduction was 
use for surface ship, and achieved good roll reduction effectiveness [6]. An adaptive backstepping controller 
and sliding mode controller were propose for active fin with radial basic function (RBF) neural network. The 
effectiveness of this device is highly dependent on the controller, therefore, researching to design an active 
fin controller with high roll reduction efficiency is also a very important task. 

Nowadays, control techniques have great progress, therefore there are many control methods have 
been researched for roll reduction control with fin stabilizer system. With advantages are reliability and 
simplicity, the proportional integral derivative (PID) control method was applied to reduce the roll motion in 
[7]. But, the disadvantage of the PID method is reduced adaptability and robustness of the system. Otherwise, 
design a controller within a linear model will be difficult to achieve high effectiveness of roll reduction due 
to roll motion is caused by the nonlinear dynamics when the ships are sailing in rough sea condition. To solve 
parameters variable and nonlinear dynamics problems, there are many robust algorithms, such as model 
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predictive control [8], [9], sliding-mode control [10], [11], fuzzy control [12], [13], were proposed for fin 
stabilizer system. However, with the ship model parameters uncertainty problem, adaptive backstepping 
control is considered as good solution for active fin system. As we all know, uncertainty parameters are 
caused by the system perturbation are compensated by using adaptive control technique, the update laws 
were used to replace the unknown parameters, and the nonlinear system mismatch problem is solved by using 
backstepping method. The virtual control input was chosen as the appropriate state variable, control variable 
output of each backstepping process was produced according to the previous process, and when ended the 
repeat process, the actual controller was obtained by optimal the Lyapunov function. 

In recent years, how to improve roll reduction efficiency received significant attention in the ship 
control field. To improve roll reduction efficiency, the main task is to overcome problems which generated 
by disturbances. With fin stabilizer system, the uncertainties mainly generated by the unknown disturbances 
which was not repeated or modeled. To reduce above mention uncertainties, there are some control methods 
which have been studied for ship active fin system such as, robust control, Hoo control [14], neural networks 
(NNs) control [15]. However, in the modern manufacture, the faster and accurate response is requested even 
under larger disturbances. To solve this problem, disturbance observer (DOB) was applied in [16], [17]. The 
approximate disturbance is estimated by comparing the real system’s control input and the symbolic model is 
main idea of the DOB. Through the symbolic model output response, the control input signal of virtual 
system is defined. The symbolic system is made by estimated disturbances, which is the feedback signal as a 
cancellation signal. DOB has some advantages can be named as robustness, less computational burden and 
high performance. For example in [18], based on DOB, a robust leader follower synchronization was 
proposed for under-actuated vessels. The DOB is effectively applied to the ship surge and yaw motions in 
this study and achieved good disturbance estimation. In our work, the DOB will be applied to the roll motion 
of the ship. In [19], an adaptive synchronization was proposed for marine surface vessels based on DOB, in 
this research, the leader-follower ship synchronization anti disturbance capability was improved by using 
DOB. Studies have also shown that disturbance observer also achieves good results with fast time varying 
disturbances, although it is designed for slow variable disturbances. 

As we all known, NNs is widely researched and achieved good results in ship control field. in [20], 
the functions parametric perturbation of the fin stabilizer system were tracked by used the NNs and achieved 
unexpected results. Ghassemi et al. [21] showed the controller which was combined between PID control and 
neural network for roll reduction control with low draught ship. In [22], an adaptive backstepping with RBF 
neural network was given for active fin without accurate model, the roll reduction effectiveness and the 
system robustness were obtained with the designed control strategy, however, the controller singularity can 
cause a potential risk. In our work, by combining the adaptive backstepping method with NNs and DOB, the 
proposed method can overcome the problems of controller singularity and mismatch, the perturbation 
parameters of the active fin system were also overcome by used the mentioned above methods. 

An adaptive controller for ship nonlinear active fin system with disturbances is presented in this 
paper. The proposed controller is combination the adaptive backstepping method with the radial basic 
function NN and the disturbance observer technique. The given scheme guarantees the stability of all closed- 
loop signals. With proposed controller, the system parameters uncertainty can be solved, the system 
disturbances can be estimated, and this method overcomes the mismatch problems, the response time is 
greatly reduced, and will achieve better roll reduction robustness and effectiveness. 

The paper’s rest are arranged as: the ship nonlinear active fin stabilizer system, wave disturbance 
model and disturbance observer will be presented in section 2. Controller design is proposed in section 3. 
Section 4 is simulation results. Finally, section 5 is this paper concludes. 


2. PRELIMINARIES 

When sailing in rough sea conditions, the ship stability is affected by roll motion, determine exact 
dynamics forces that affect to the vessel motions is very difficult, in this paper, we study the stability of the 
vessel in relative to roll motion. So, the ship nonlinear active fin model is given as [23], [24]. 


(Mx t AM,)@ a gnP aa FwHlPl T Dphmo[1 Eä (~/F51)*] = Het + Awe (1) 
In (1), 
My : inertia roll moment; 
AM; : added inertia roll moment; 
p : ship’s roll angle; 
p : ship’s roll rate; 
Pri : ship’s flooded angle; 
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Hot : fin’s control moment; 

Hwt : sea wave’s acting moment; 

Dy : displacement of the ship; 

hm : ship’s initial metacentric height; 


čn and čų : damping coefficients. 
Mx, AM, n» Ey, Het can be given as: 


pee 0.0227B,, 0.0043L 
M, + AM, = z (0.3085 pe eee ’y A 
rf 


_ 2rıyDphm(M, + 4M) , _3r,(Mx + 4M,) 
n7 T ow = 4 (3) 


gly 
Het = —pv* Spl Cis (Gp + ee (4) 


where B,e is the ship’s breadth, Lp is the ship’s length, d,- is the vessel draft, r1, r) are test factors, pis 
seawater density, v is the vessel speed, Sp is fin area, ly is the fin’s acting force arm, Cj, is the lift coefficient 
of fin, and aș is the rotation angle of fin. Take AM,, My, Gy, & and Het instead of (1), we can get: 


Ë =P + 729° +Y% + YaGlG| + oas +d (5) 


where y1, Y2, Y3, Y4, O are coefficients, d is sea wave disturbance. 
Denoted the state variables x = [x,,x2]" = [p,@]", the output signal y = g = x, , and u = aç is 
the input signal. From (5) can be written as: 


X= X3 
X= f(x) + g(x)ut+d 
y =X (6) 


Assumption 1. gis an unknown smooth function, and signs of g is known. And exist constant g, > gz > 0 
satisfy g4 > |g| > gz, and Exist constant gg > 0 such that |ġil < Jia- 

Lemma 1. If V(t, x) is positive definite, and V < —k,V + kz, with kı > 0, k, > 0 are bounded, and then 
V(t, x) < kz/kı + (V (0) — kz /k;1)e "t 


2.1. Disturbance observer 
The second function in (6) is rewritten as: 


d =X, — f (x) - g&u (7) 
The disturbance observer is proposed as: 


d = -K (x1, x2)Â + K (x1, %2)(% — fŒ) - g@u) (8) 


where K (x1, xz) = a is the parameter with a > 0 is designed. Due to the disturbance’s derivative there is no 
prior information, so assuming that disturbance in (6) varies slowly [25]. 


d=0 (9) 
We defined observer error as difference between estimated and actual disturbance 


d=d-d (10) 


which can be written as in (11). 


d = —K(x,,x,)d (11) 
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We can see that, in (6), x2 is hard to obtain, this signal is essential to realize the disturbance 
observer. To solve this problem, an auxiliary variable is defined as: 


r = d—4q(x2). Defining 

d (x2) = K (x1, 2) Xo (12) 
then 

t = =K (x1, x2)r + K(x, %2)(—f Œ) — gu — q(x2)) (13) 


The estimated disturbance is obtained by: 


â =r + q(x) (14) 


and it is known that 


d=r+q(x2) (15) 
We design 
Ug = d/g(%) (16) 


The second function in (6) is rewritten as: 


i = f@)+g@)v—u)td=fX)+g@vi+d (17) 


where d = d — À, u = v — ug and the system (6) is transformed as: 


ži =x% 
4 =f) +g@)vt+d 
y=xX (18) 


The disturbance observer’s structure is shown in Figure 1. In the next section, a new design procedure is 
developed for the system (18) by using the adaptive backstepping method and radial basis function NN. 


Disturbance 
observer 
u 


Figure 1. Disturbance observer’s structure 


3. CONTROL SYSTEM DESIGN 
3.1. RBF neural network 

In nearly years, radial basic function neural network (RBFNN) have been widely studied, please 
refer [26] for detailed treatment of RBFNN. In this paper, RBFNN model is given as: 


mm (Z) = W" F(z) (19) 
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In (19), input vector z E N ER™, W = [w1, w3, ..., Wn] is weight vector, n is node number of the NN, and 
F(z) = [f, (2), (©), .--, fa (Z)] with f;(z) is chosen as Gaussian functions 


= -E-D Z-Y) 
fila) = exp a] (20) 


where i = 1,2,...,n, Yi = vin, Vias Yim] is the receptive field center, £; is the Gaussian function’s width. 
As we all known that, using RBFNN to approximate an unknown continuous function f(x) over 
compact sets, we can obtain 


fx) =W*F(x)+n Wx EQXQSR™ (21) 


In (21) the approximation error isn, the basic function F(x), and w* is ideal weight matrix. The ideal matrix 
w* is an “artificial” quantity which is given for analytical purpose and have form as: 


w* = arg min, {suplf (x) -W'F œI} (22) 
3.2. Control design 

Step 1: The variable error is defined z, = x,, the virtual control law is choosen as az. The second 
variable error is defined 

Z2 = X2 — Q2 (23) 
The z, first derivative is calculated as: 

żŻ =Z,+ az (24) 
Virtual control law is defined as: 

Q2 = CZ (25) 
The Lyapunov function is considered V, = z? /2 and the first derivative of V} as: 

V, = Ge? + Z122 (26) 

Step 2: We define Zz} = x2 — æ, and the differentiate z} with time yield’s respect is 

ż, = f(#) + gv + d — å, (27) 


with f(x) and g(x) are unknown, &, is unknown scalar function. We define h(x) and using a radial basis 
function neural network Y*T p(x) to approximate h(x) as: 


WR) = TS (F@ +d - á) = wpe) + 5° (28) 
In (28) "is the ideal weights, and the approximate error ¢* > Ois a constant, then (27) can be writen as: 
Zz = JE + YPE) +57) (29) 
The control law is chosen as: 
v = —C2Z, — W" p(X) (30) 


In (30) Wis the estimate parameter of y* matrix, denote p = W — y* is the parameter error, we have: 


Z = g(%)(—C222 — W" p(X) + ç*) (31) 
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Vz Sy 
1" 29%) 


+50'T 1p 


with T = T7 > 0 is the adaptation gain vector, the first derivetive of V, is 


, , 5 en EE: 23.g(%) 
V, = V, — 0225 + 226° + WT" ($ — Tpz2) + 2(g(%))? 


The adaptation law is considered as in (34): 


$ = T (pz — op) 


o 


(32) 


(33) 


(34) 


with ô > 0 is a constant. By using (26), (31), and (34), we obtain the first derivative of V, as in (35). 


; ann 25gG(X) 
V, = c,22 — C923 + 2,2, + 26° — Op p+ CIEE 


Refer to the facts 


1 1 42 
ZZ < Z? +223; 26° < zz2s" 


then 
5 1a 
-(& -1 z? - (c> -2)2 +3 i 
(cy Zz 4) 72 4° 
z3G(%) 
-2 - + 72d 
2s II’ lp; |? ) 2(g(x))? 
Note that, 
"ORR 
ee ae Ce ee 4 Aga) 
(223 + Faq) E (6222 + age) 


Cz is choosen such that 


« _ (C2-98) 
2 Rw 7 


Denote that 
Cy = 1 = c1; CQ — : = Caig” +Š yill? =k 

We pay attention to |ç*| < çy and ||Y*|| < Yy. Then can get: 
kisi + yh = ky = L 


then 


V, < — 2A 2 a 
2 5 — C1171 — C2222 eat 
1 


Let 


Cc ô 
z; C 
29(x) 24-1max5 


Cc. aan 
1022 = 


C. — 
11 2 
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(37) 


(38) 


(39) 


(40) 


(41) 


(42) 


(43) 
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with C > 0 is bound constant and satisfies 6 > CA — 1max and it follows 


V <—CV, +L (44) 


with the assumptions and lemma 1 which are given in section 2, (44) shows that V, is bounded. That means 
the bounded stability of system (6) is guaranteed. 


4. SIMULATION RESEARCH 

The design vessel is used for simulation, which have the parameters as given in Table 1. In the 
simulation, the parameters were selected as: yı = —0.2117, yz = 0.2939, y3 = 0.6748, y, = —0.113, 
ô = 0.01, T = diag{0.001}, c = 0.1, c, = 2, fin angle |a| < 25deg. n = 135 is node number of the NN, 
€; = 2(i = 1...25). External disturbance is given as square-wave and sinusoidal signal. The control 
performance is illustrated in below simulation results. The disturbance is choice as sinusoidal signal d = 
0.25 sin(0.8t) and square-wave signal: d = 0.25(sign(sin(0.8t))). 


Table 1. Vessel parameters 
Vessel parameters 


Ship’s length 175 m 
Ship’s breadth 25m 
Draft 8.0 m 
Displacement 21130 t 
Fin area 20 m? 
The fin acting force arm 15.0 m 
Fin fift coefficient 3.4 
Ship’s flooding angle 45 (°) 
Shp’s Initial metacentric height 5.5m 
Ship’s speed 18 kts 


As we shown in Figure 2 and Figure 3 the roll angle is reduced to the desire value and this value is 
maintained that mean the effectiveness of the proposed controller is illustrated. Figure 4 and Figure 5 shown 
that the system’s disturbance is estimated very well with proposed disturbance observer. Figure 6 and 
Figure 7 shown that the controller designed has the response time close to 2 s while response time of the 
controller without disturbance observer close to 10s which is shown in Figure 8 and close to 20 s which is 
shown in Figure 9, that mean, by using disturbance observer, the response time is greatly reduced. 


25 T T T T T 25 T T T T T 
—— — Tracking angle/degree ——— Tracking angle 
Roll angle/degree 20 Roll angle 

g 3 
4 4 
© T 

5 L 4 l 4 l 15 f i L f f 

0 5 10 15 20 25 30 5 10 15 0 5 30 
Time/s i 
Time/s 
Figure 2. The ship roll angle with Figure 3. The ship roll angle with 
d = 0.25 sin( 0.8 t) d = 0.25(sign(sin( 0.8 t))) 
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© 
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-0.1 -0.44 
-0.2 06} 
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Figure 4. Disturbance observer with Figure 5. Disturbance observer with 
d = 0.25 sin( 0.8 t) d = 0.25(sign(sin( 0.8 t))) 
30 : : : r 30 
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20+ J 
20H 
15} 
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-20 1 L 1 L 1 -30 L 4 L 1 1 
a 5 10 15 20 25 30 a 5 10 15 20 25 30 
Time/s Time/s 
Figure 6. The fin angle with Figure 7. The fin angle with 
d = 0.25 sin( 0.8 t) d = 0.25(sign(sin( 0.8 t))) 
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Fin angle without disturbance observer Fin angle without disturbance observer 
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20 
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Figure 8. The fin angle without disturbance observer 
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with d = 0.25 sin( 0.8 t) 


L 
0 10 20 30 40 50 60 70 8 90 
Time/s 


Figure 9. The fin angle without disturbance 
observer with d = 0.25(sign(sin( 0.8 t))) 
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5. CONCLUSION 

In our work, an adaptive control scheme based on NN and disturbance observer for ship nonlinear 
active fin system is proposed. One disturbance observer is given to estimate the system’s disturbances, in 
addition, RBFNN is proposed to approach the unknown elements in the ship nonlinear active fin system, and 
therefore the system can obtain good roll reduction effectiveness and overcome the uncertainties of the 
model. Unknown disturbance and system parameters were solved with the proposed method. The simulation 
results illustrate that good roll reduction efficiency can be achieved. 
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